N88-12754 


NASA CONTRACTOR REPORT 

(NASA-CE- 179222) A MOITII1E-U RE-SCALE 
7CBBDIEKCE MODEL EASED CN VAFIAEIE 
FAHTITICNING CF 1 CRBOIFHT KISITIC ENEBGI 
EFFCTBUM Final Report (Uoivc rsities Space Unclas 

Fesearcb Association) 42 F CSCL 20D G3/34 0111209 

NASA CR-179222 


A MULTIPLE-TIME-SCALE TURBULENCE MODEL BASED ON VARIABLE 
PARTITIONING OF TURBULENT KINETIC ENERGY SPECTRUM 


By S. -W. Kim 

Universities Space Research Association 
Systems Dynamics Laboratory 
Science and Engineering Directorate 

and 


C. -P. Chen 

Department of Mechanical Engineering 
University of Alabama in Huntsville 
Huntsville, Alabama 35899 


Final Report 


November 1987 


Prepared for 

NASA- Mar shall Space Flight Center 

Marshall Space Flight Center, Alabama 35812 


r 


L 


l* 




TECHNICAL REPORT STANDARD TITLE PAGE 


1. REPORT NO. 2. GOVERNMENT ACCESSION NO. 

NASA CR- 179222 


4 . TITLE AND SUBTITLE 

A Multiple- Time- Scale Turbulence Model Based on Variable 
Partitioning of Turbulent Kinetic Energy Spectrum 


7. AUTHOR(S) 

S.-W. Kim* and C.-P. Chen** 


9. PERFORMING ORGANIZATION NAME AND ADDRESS 


3. RECIPIENT'S CATALOG NO. 


5. REPORT DATE 

November 1987 


6. PERFORMING ORGANIZATION CODE 


6. PERFORMING ORGANIZATION REPORT » 


|10. WORK UNIT NO. 


George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 


12. SPONSORING AGENCY NAME ANO ADORESS 


National Aeronautics and Space Administration 
Washington, D.C. 20546 


15. SUPPLEMENTARY NOTES 


11. CONTRACT OR GRANT NO. 

NAS8- 35918 


13. TYPE OF REPOR'i & PERIOD COVERED 

Contractor Report 


14. SPONSORING AGENCY CODE 


Contract Monitor: N. C. Costes , Systems Dynamics Lab., 

Marshall Space Flight Center 
ice Research Association 


*Universities Space Research Association 
**Dept. of Mechanical Engineering, University of Alabama in Huntsville. 


16. ABSTRACT 


A multiple- time- scale turbulence model of a single point closure and a simplified 
split- spectrum method is presented in this paper. In the model, the effect of the 
ratio of the production rate to the dissipation rate on eddy viscosity is modeled by 
use of the multiple -time -scales and a variable partitioning of the turbulent kinetic 
energy spectrum . The concept of a variable partitioning of the turbulent kinetic 
energy spectrum and the rest of the model details are based on the previously 
reported algebraic stress turbulence model. Example problems considered include: 
a fully developed channel flow, a plane jet exhausting into a moving stream, a wall 
jet flow, and a weekly coupled wake-boundary layer interaction flow. The computa- 
tional results compared favorably with those obtained by using the algebraic stress 
turbulence model as well as experimental data. The present turbulence model, as 
well as the algebraic stress turbulence model, yielded significantly improved compu- 
tational results for the complex turbulent boundary layer flows, such as the wall 
jet flow and the wake-boundary layer interaction flow, compared with available 
computational results obtained by using the standard k- e turbulence model. 


17. keywords Multiple- Time- Scale Turbulence ,8 - distribution statement 
Model, Algebraic Stress Turbulence Model, 

Turbulent Wall Jet, Turbulent Wake- Unclassifi 

Boundary Layer Interaction Flow , Turbulent 

Backward- Facing Step Flow, Confined 

Swirling Jet, Finite Element Method, Control- 

Volume Based Finite Difference Method 


19. SECURITY CLASSIF. (of thi. r.port> 

Unclassified 


Unclassified /Unlimited 


MS PC - Form 3192 (May 1969) 


■ 

20. SECURITY CLASSIF. (of thi. page) 

21. NO. OF PAGES 

22. PRICE 

Unclassified 

41 

NTIS 

! 




For tale by National Technical Information Service, Springfield, Virginia 221 51 














TABLE OF CONTENTS 


I. INTRODUCTION 

II. MULTIPLE-TIME-SCALE TURBULENCE MODEL 

2-1. Turbulence Equations 

2-2. Eddy Viscosity Equation 

2- 3. Wall Function Boundary Conditions 

III. COMPUTATIONAL RESULTS 

3- 1. A Fully Developed Channel Flow 

3-2. A Plane Jet Exhausting into a Moving Stream 
3-3. A Wall Jet Issuing into a Moving Stream — ., 
3-4. A Wake-Boundary Layer Interaction Flow .... 

IV. CONCLUSIONS AND DISCUSSION 


APPENDIX I. THE TURBULENCE MODEL CONSTANTS 


A- 1-1. Free Stream Turbulence 

A- 1-2. Near Wall Equilibrium Flow 

A- 1-3. Homogeneous Turbulent Shear Flow 


APPENDIX II. A CONTROL- VOLUME BASED FINITE DIFFERENCE 
COMPUTATION OF SEPARATED AND SWIRLING 
TURBULENT FLOWS 


A-2-1. A Turbulent Backward -Facing Step Flow 
A- 2- 2. A Confined Swirling Jet 


REFERENCES 


Page 

1 

2 

2 

4 

5 

6 

6 

8 

11 

14 

16 

17 

17 

19 

20 


23 

23 

28 

34 


PRECEDING PAGE BLANK 


NOT KILLED 


iii 


LIST OF ILLUSTRATIONS 


Figure Title Page 

1. Description and nomenclature of the multiple-time- scale 

turbulence model 2 

2. A fully developed channel flow 7 


(a) Velocity, turbulent kinetic energy, and Reynolds stress 

(b) Dissipation rate, k t /k p , and e t /e 

3. A plane jet exhausting into a moving stream 8 

(a) Development of the jet along the flow direction 

(b) Velocity, turbulent kinetic energy, and Reynolds stress 
at x=0.381 m 

(c) Dissipation rate, k./k , and e./£ at x=0.381 m 

(d) Velocity, turbulent kinetic energy and Reynolds stress 
at x=0.650 m 

(e) Dissipation rate, k t /k p , and e t /e at x=0.650 m 

4. A wall- jet flow 11 

(a) Development of the wall-jet flow along the flow direction 

(b) Velocity, turbulent kinetic energy, and Reynolds stress 
at x=l. 1532 m 

(c) Dissipation rate, k./k , and e./e at x=1.1532 m 

t p l P 

(d) Velocity, turbulent kinetic energy, and Reynolds stress 
at x=1.6882 m 

(e) Dissipation rate, k./k and e./e at x=1.6882 m 

t p t p 

5. A wake-boundary layer interaction flow., 14 

(a) Velocity, turbulent kinetic energy, and Reynolds stress 
at x=0.40 m 

(b) Dissipation rate, k /k , and e./e at x=0.40 m 

x p x p 

(c) Velocity, turbulent kinetic energy, and Reynolds stress 


at x=0.86 m 

(d) Dissipation rate, k./k , and e + /e_ at x=0.86 m 

x p x p 

A . 1 Streamline contour for backward- facing step flow 24 

A. 2. Turbulent kinetic energy contour for backward- facing 

step flow 24 

A. 3. Velocity profiles for backward- facing step flow 25 

A. 4. Turbulent kinetic energy profiles for backward- facing 

step flow 25 

A. 5. Reynolds stress profiles for backward- facing step flow 26 


IV 



LIST OF ILLUSTRATIONS (Concluded) 


Figure Title 

A. 6. The ratio of dissipation rate to energy transfer rate for the 
backward- facing step flow 


A. 7. 
A. 8. 
A. 9. 


Streamline contour for backward- facing step flow, 

o. = a, . = 1 . 0 and a = a , = 1.3 

kp kt ep et 

Turbulent kinetic energy contour for backward- facing step flow, 

o. =<*..=1.0 and <j = a . = 1.3 

kp kt ep et 

Streamline contour for confined swirling jet 


A. 10. Turbulent kinetic energy contour for confined swirling jet 

A. 11. Development of the axial mean velocity along the center line for 
confined swirling jet 

A. 12. Axial velocity profiles of the confined swirling jet 


A. 13. Swirl velocity profiles of the confined swirling jet 

A. 14. R.M.S. values of the fluctuating axial velocity 

A. 15. The ratio of dissipation rate to energy transfer rate for 
confined swirling jet 


A. 16. 
A. 17. 


Streamline contour for confined swirling jet, 

o. = cr. , = 1.0 and a _ - a, = 1.3. ....... 

kp kt ep et 


Turbulent 
°kp = a kt 


kinetic energy contour for confined swirling jet , 

= 1.0 and a = o. = 1.3 

ep et 


Page 

26 

27 

27 

29 

29 

29 

30 

30 

31 

31 

32 
32 


v 


LIST OF TABLES 


Table Title Page 

A.l. Contour Labels for Backward- Facing Step Flow 33 

(a) Streamline (ip) 

(b) Turbulent Kinetic Energy (k) 

A. 2. Contour Labels for Confined Swirling Jet 33 

(a) Streamline (\p) 

(b) Turbulent Kinetic Energy (k) 



CONTRACTOR REPORT 


A MULTIPLE -TIME -SCALE TURBULENCE MODEL BASED ON VARIABLE 
PARTITIONING OF TURBULENT KINETIC ENERGY SPECTRUM 

I. INTRODUCTION 


In turbulent flows, instability of mean flow generates large eddies, the large 
eddies break up into smaller eddies , and the fine scale eddies are dissipated by 
viscous force. The multiple- time- scale turbulence model is derived by partitioning 
the turbulent kinetic energy spectrum into a set of wave number regions, integrating 
the partitioned energy spectral density equation to obtain energy evolution for each 
of the wave number regions, and considering the energy flux rate between adjacent 
wave number regions. Therefore, the generation, cascade, and dissipation of turbu- 
lent kinetic energy are considered in the multiple-time-scale turbulence model; 
whereas, only the generation and the dissipation of turbulent kinetie energy are 
considered in the classical turbulence models such as the k- e turbulence models and 
the Reynolds stress turbulence models. 

A few papers on the multiple-time- scale turbulence model have appeared 
recently. The multiple-time- scale turbulence model used in Hanjelic et al. [1] and 
Fabris et al. [2] are based on the simplified split- spectrum method. Generalization 
of the model for multiple split- spectrum case and extension of the multiple-time-scale 
concept to convection- diffusion of scalar variables is given in Schiestel [3,4] . The 
multiple-time-scale turbulence model presented herein is based on the simplified split- 
spectrum method and variable partition of the turbulent kinetic spectrum. The rest 
of the model details are based on the algebraic stress turbulence model given in 
References 5 through 7, among many others. 

Uniform division of the turbulent kinetic energy spectrum has been used in 
the multiple-time- scale turbulence models proposed in References 1 through 4. A 
difficulty may arise in this case, since the turbulence field of complex turbulent flows 
is characterized by inhomogeneous turbulent kinetic energy spectrum distributed in a 
wide range of wave numbers. If the partition is moved to too high wave numbers, 
then the multiple- scale turbulence model may reduce to a single scale model, on the 
other hand, if the partition is moved to too low wave numbers, then production of 
turbulent kinetic energy will be contained in the dissipation range which is contra- 
dictory to the multiple-time-scale concept [1-4]. This difficulty can be eliminated by 
using a variable partitioning of the turbulent kinetic energy spectrum in such a way 
that the partition is moved to the high wave number direction when production is 
high, and that the partition is moved to the low wave number direction when produc- 
tion vanishes. In the present turbulence mdoel, the partitioning of the turbulent 
kinetic energy spectrum is dependent on the turbulence intensities, production rate, 
energy transfer rate:-, and dissipation rate, and is determined as a part of the solu- 
tion. Furthermore, the variable partitioning method rendered the effective eddy 
viscosity coefficient to decrease when produciton is high and to increase when produc- 
tion vanishes. The same effect of the production rate on the eddy viscosity coeffi- 
cient can be found in experimental data [8] as well as in algebraic stress turbulence 
models [5-7], 

The turbulent boundary layer flows considered herein were solved by a finite 
element method [5,9,10], It has been shown in Reference 10 that the finite element 


method could solve a wide range of laminar boundary layer flows, such as the Blasius 
flat plate flow, the retarded Howarth flow, flow over a wedge, plane stagnation flow, 
flow over a circular cylinder, flow in the wake of a flat plate, uniform suction flow 
over a flat plate, flow over a cone, and flow over a sphere, as accurately as any 
available numerical methods including the semi- analytical methods. The numerical 
details on turbulence equations can be found in Reference 5. 

Establishment of the turbulence model constants is discussed in Appendix I. 

Use of the present turbulence model in finite difference computation of elliptic turbu- 
lent flows improved the computational results in comparison with those obtained by 
using the standard k-e turbulence model as can be found in Appendix II. 


II. MULTIPLE-TIME-SCALE TURBULENCE MODEL 


The turbulent boundary layer flow equations are given as: 
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where u and v are the time averaged velocities in flow direction and in transverse 
direction, respectively, p is the pressure, v is the kinematic viscosity of the fluid, 
and is the turbulent eddy viscosity. 


2-1. Turbulence Equations 

The underlying idea [1-4] and the nomenclature of the multiple- scale turbulence 
model is shown in Figure 1, where k^ is the turbulent kinetic energy of eddies in the 

production range, k^. is the turbulent kinetic energy of eddies in the dissipation range, 
Pr is the production rate of turbulent kinetic energy, k (k=kp+kp is the total turbu- 
lent kinetic energy, e is the transfer rate of turbulent kinetic energy from the pro- 

r 

duction range to the dissipation range, and is the dissipation rate. 
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Figure 1. Description and nomenclature of the multiple-time-scale turbulence model. 
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Transport equations for the turbulent kinetic energies are given as [1 4]: 


u 




9 x 


ay 


’y v 


„ ] = Pr - e 

9y p 


( 3 ) 


3k t sk t a „ v t . 3k t. 

Tx" v Ty - ’ Jy 1(v + ^ TjT 


= £ 


(4) 


where a, and a, . are constant coefficients, 
kp kt 

The convection- diffusion equations for the energy transfer rate and the dissipa- 
tion rate can be written as, in general: 
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Following the work of Hanjelic et al. [1] and based on the algebraic stress tur 
bulence model [5-6], the influence functions proposed herein are given as f = 

C pl Pr2/k P + C p2 Pl V k p ' C p3 £ p 2/k P and f t = C tl E p 2,k t + C t2 e p^t /k t ' C t3 e t 2/k f 
respectively. The complete form of the convection- diffusion equations for the energy 
transfer rate and the dissipation rate are given as: 
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where a and are constant coefficients, c pJ ^ (£-1,3) and c^ (£-1,3) are turbu- 
lence model constants. The c pl Pr 2 /k p and e u e p 2 /k t terms were obtained from 

physical dimensional analysis. The first term increases the energy transfer rate when 
production is high , and the second term increases the dissipation rate when the energy 
transfer rate is high. Both of these terms are equivalent to the variable energy 
transfer function proposed in Hanjelic et al. [1], which is given as -0. 3(k p /k^.- 1) / 


( V k t + 1 V /k P ' 


It can be shown that the variable energy transfer function proposed 


3 


in Hanjalic et al. is equivalent to those of the present turbulence model, if the ratio of 
kp to k t approaches a constant which is determined by the ratio of production rate to 

dissipation rate as discussed in Appendix I. Further discussion on the variable 
energy transfer function can be found in Pope [11] and Hanjalic and Launder [12], 
among many others. 


°kt = 
and c 


The model constants used in the present study are; 
0.75, a £t = 1.15, c pl = 0.21, c p2 = 1.24, c p3 = 1.84, 
t3 = 1 . 66 . 



0.75, a = 1.15, 

ep 

0.29, c t2 = 1.28, 


2. Eddy Vis*. si+' Equation 

The eddy viscosity expression proposed herein is given as: 


v 


t 



(9) 


where c^ (=0.09) is a constant. The eddy viscosity equation is equivalent to that 

of the standard k-e turbulence model except that the dissipation rate (e) has been 
replaced by the energy transfer rate of large eddies. Therefore, the turbulence 
length scale is related to the energy transfer rate (e ), rather than the dissipation 
rate (e^. P 

Experimental data shows that the eddy viscosity coefficient (c ) decreases as 

the ratio of Pr/e^ is increased; and increases as the ratio of Pr/e is decreased. In 

algebraic stress turbulence models, the eddy viscosity coefficient (c ) is given as a 

function of the ratio of production rate to dissipation rate [5-7]. In the present 
multiple-time-scale turbulence model, which is based on the variable partitioning of 
the turbulent kinetic energy spectrum method, the c function is implicitly included 
in the eddy viscosity equation, i.e., y 


v t ' c „f k2/E P * ° 1J k2/E t <10 > 


where = c 1J f G ^/ £ p is defined to be the effective eddy viscosity coefficient in the 
following discussions. The ratio of decreases as the ratio of Pr/e^. is increased; 

and increases as the ratio of Pr/e^. is decreased. In the near wall regions, the pro- 
duction rate is approximately equal to the dissipation rate, thus the energy transfer 
rate has to be approximately equal to both of these. In the free stream region where 
production vanishes, the turbulence model constants given previously yield the value 
°f e t^ e p to be 2 * 5 approximately (see Appendix I). 

The dissipation rate is weakly coupled with the rest of the turbulence quantities 
in the multiple-time-scale turbulence model than in the single -time- scale turbulence 
models such as the k-e type turbulence models. If the multiple -time- scale concept is 
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to be justified, then it is expected that the computational result of the dissipation 
rate (e^) obtained by using the multiple- time- scale turbulence model would be com- 
parable to that obtained by using the single -time- scale turbulence model such as the 
algebraic stress turbulence model [5-6]. The present computational results showed 
that the above requirement was very well satisfied. 

The eddy viscosity equation proposed by Hanjalic et al., which is given as: 
v . — c (k +kjk /e , is compatible with the near wall mixing length theory or the 

t y P t p p t 

standard wall function method only when k t vanishes in the near wall region. On the 

other hand, the eddy viscosity expression given in equation (9) is compatible with 
both the near wall mixing length theory and the wall function method in general. 


2-3. Wall Function Boundary Conditions 

The wall function boundary conditions were derived from the standard wall 
function method [13]. These are given as: 


u 1 . + N 

— = - In (Ey ) 
u k J 


( 11 ) 


x = - [ p kc 1/4 k 1/2 /ln(Ey + ) ] u 
w y 


k = 


c 

y 



/ p 


e 


1 


c 3/4 k 3/2 

y 



K 

1/2 ; 

CI c (c «, 
ep V p3 


c . - c 0 ) 
pi p2 


( 12 ) 

(13) 

(14) 

(15) 



(16) 


where u + , u + =u/u , is a non-dimensional velocity, u , u = /(x / p), is the wall fric- 

T XT W 

+ + . 

tion velocity; x is the wall shearing stress; y , y = u y/v, is the wall coordinate; 

k is the von Karmann constant; and E is an experimentally determined constant coef- 
ficient. k= 0.41 and E=9.0 has been used in the present study. Derivation of equa- 
tion ( 15) is given in Appendix I . 
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III. COMPUTATIONAL RESULTS 


The governing differential equations were solved on the physical domain using 

3 

physical dimensions. Whenever necessary, the values of 1.225 kg/m and 0.17854x 
10 ^ kg/m- sec were used for density and molecular viscosity, respectively. 

The initial condition data for velocity, turbulent kinetic energy, and dissipation 
rate were obtained from experimental data [5-6]. The initial condition data for the 
ratios of k^/k^ and £p/^ were obtained by interpolating the near wall values and the 

free stream values using an equation of the form: 


r = r e + (r w " r e ) /(Pr/e t } 


(19) 


where r stands for k /k. or e /e, , the subscript e denotes the external free stream 

P t P l 


value, the subscript w denotes the near wall value, and the producti on ra te was 
obtained from the experimental data using the relationship that Pr = -u'v' (3u/3y). 


The wall function boundary conditions, equations (11-16), were used in the near 
wall region and the vanishing gradient boundary condition was used at the outer edge 
of the computational domain for all the cases, unless otherwise specified. Details on 
the computational procedure can be found in Reference 5. The convergence criterion 
used is given as: 



n+1 




n+1 


< e 


j = 1 , N , no sum on j 


( 20 ) 


where a. stands for the nodal values of u , k , e 


) 


-7 


or N denotes the total number 


of degrees of freedom; and e = 1x10 was used. 


3-1. Fully Developed Channel Flow 

The experimental data for the fully developed channel flow considered herein 
can be found in Laufer [14]. The Reynolds number based on the channel half width 
of 0.0635 m and the center line mean velocity of 7.07 m/sec is approximately 30,800. 

The computational domain extending from y = 0.005 m near the wall, which 
corresponds to y + = 100, to y = 0.0635 m at the center of the channel was discretized 
using 20 equally spaced quadratic elements. The Dirichlet boundary conditions for u, 
kp and k^. at the near wall region were obtained from experimental data and by using 

equation ( 15) , whereas the boundary condition for the dissipation rate was obtained 
from the near wall mixing length theory [5-6]. The near wall boundary conditions 

used are: u = 5.084 m/sec; k = 0.213 m^/sec^; k = 0.172 m^/sec^; k, = 0.213 m^/ 

2 2 3 p 

sec , and = 10.64 m /sec . 
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The flat profiles for the velocity, both of the turbulent kinetic energies, the 
energy transfer rate, and the dissipation rate were used as initial guess, and the 
convergence criterion stated earlier was achieved in 150 iterations. The computational 
results for the velocity, the turbulent kinetic energy, and the Reynolds stress 
obtained by using the multiple-time-scale turbulence model are compared with those 
obtained by using the algebraic stress turbulence model as well as experimental data 
in Figure 2(a). It can be seen that the multiple-time- scale turbulence model yielded 
almost the same computational results as the algebraic stress turbulence model. The 
dissipation rate obtained by using the multiple- time- scale turbulence model was found 
to be the same as that of the algebraic stress turbulence model. The ratios of k t /k 


and e t /e are shown in Figure 2(b). It can be seen in Figure 2(b) that the value of 
e^/e is approximately equal to unity, and hence the effective eddy viscosity coeffi- 
cient (c^ ) is approximately equal to 0.09 in the entire flow domain. 
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ENERGY, AND REYNOLDS STRESS, x. o, A : EXPT DATA FOR VELOCITY, TURBULENT KINETIC 
ENERGY, ANO REYNOLDS STRESS. M: MULTIPLE-TIME-SCALE TURBULENCE MOOEL, 

A: ALGEBRAIC STRESS TURBULENCE MODEL 
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Figure 2. A fully developed channel flow, (a) Velocity, turbulent kinetic energy, 
and Reynolds stress, (b) Dissipation rate, k^/kp an d e^/Cp. 
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3-2. A Plane Jet Exhausting into a Moving Stream 

A plane jet exhausting into a moving stream is considered below, the experi- 
mental data of which can be found in Bradbury [15]. In the following discussions, 
the half jet width is defined as a distance from the center line of jet to a location 

where the excess velocity is half of the center line excess velocity. 

The computational domain in flow direction starting from x = 0.095 m (x/d = 10, 
d = 0.009525 m is the jet exit width) to x = 0.65 m (x/d = 70) was discretized by 

580 line-steps; and the transverse domain extending y = 0 m (center line of the jet) 

to y = 0.1 m (approximately 12.5 times of the half jet width at x/d = 10) was discre- 
tized by 45 unequally spaced quadratic elements. 

Decay of the center line velocity and evolution of the half jet width along the 
flow direction, obtained by using the two turbulence models, compared favorably 
with experimental data as shown in Figure 3(a). The computed velocity profile, 
turbulent kinetic energy profile, and the Reynolds stress are compared with experi- 
mental data in Figures 3(b) and 3(d) at two flow -direction locations. These profiles 
compare more favorably with experimental data at the middle of the flow direction 
domain than at the downstream end of the computational domain. It can be seen that 
the multiple- time -scale turbulence model yielded slightly better computational results 
than the algebraic stress turbulence model at the downstream end of the computational 
domain. The dissipation rate obtained by using the multiple- time- scale turbulence 
model compared favorably with that of the algebraic stress turbulence model as shown 
in Figures 3(c) and 3(e). The ratio of e^/e exhibited sharp transition in the viscous 

super layer (the transition region between the high turbulence region and the free 
stream low turbulence region), thus the effective eddy viscosity coefficient (c^) also 

changed abruptly, but still smoothly and continuously, in the viscous super layer. 

An average of 9 iterations was required to satisfy the previously stated con- 
vergence criterion for each line- step. 


y 1/2 uo 
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(a) Development of the jet along the flow direction. 
Figure 3. A plane jet exhausting into a moving stream. 
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ENERGY, AND REYNOLDS STRESS, M: MULTIPLE-TIME-SCALE TURBULENCE MODEL, 

A: ALGEBRAIC STRESS TURBULENCE MODEL 

(b) Velocity, turbulent kinetic energy and Reynolds stress at x = 0,381 m. 
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(c) Dissipation rate, k t /kp, and at x = 0.381 m. 


Figure 3. (Continued) 
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, : COMPUT. RESULTS FOR VELOCITY. TURBULENT KINETIC 

ENERGY, AND REYNOLDS STRESS, x, o.A : EXP'T DATA FOR VELOCITY, TURBULENT KINETIC 
ENERGY, AND REYNOLDS STRESS, M: MULTIPLE-TIME-SCALE TURBULENCE MODEL, 

A: ALGEBRAIC STRESS TURBULENCE MODEL. 

(d) Velocity, turbulent kinetic energy, and Reynolds stress at x = 0.650 m. 
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(e) Dissipation rate, k^/kp, and e^/e at x = 0.650 m. 
Figure 3. (Concluded). 
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3-3. A Wall Jet Issuing into a Moving Stream 

One of the most complete experimental data for the wall jet flows can be found 
in Irwin [ 16] , and a finite difference computation of wall jet flows can be found in 
Ljuboja and Rodi [17], among many others. 

Input data used in computation of the flow were obtained directly and/or by 
curve- fitting the experimental data [5-6]. The computational domain extended from 
y = 0.002 m to y = 0.08 m in the transverse direction; and from x = 0.5532 m (x/b = 
82.2) to x = 1.6892 m (x/b = 251) in the flow direction. The transverse domain was 
discretized by 45 unequally spaced quadratic elements, and the flow direction domain 
was discretized by 1135 line- steps. An average of 12 iterations were required for 
each line-step to achieve the same convergence criterion given previously. 

The computational results of the flow development along the downstream direc- 
tion obtained by using the multiple-time- scale turbulence model are compared with 
that of the algebraic stress turbulence model as well as experimental data in Figure 
4(a). It was found that the present turbulence model yielded significantly improved 
computational results compared with those of Reference 17. The velocity, turbulent 
kinetic energy, and Reynolds stress at the two downstream locations are shown in 
Figures 4(b) and 4(d). The dissipation rates obtained by using the two turbulence 
models were found to be almost identical for the wall jet flow ; and the ratio of e t / e p 

exhibited sharp transition in the viscous super layer [see Figures 4(c) and 4(e)]. 
Both of the turbulence models yielded almost identical computational results for the 
wall- jet flow; nevertheless, the multiple -time -scale turbulence model yielded slightly 
better computational results in the outer edge of the boundary layer. 

y 1/2 T w u o u e 



, , : COMPUT MAX. EXCESS VELOCITY (u 0 ). EXTERNAL VELOCITY (u e ), 

WALL SHEARING STRESS (t*). AND LENGTH SCALE (y 1/2 ), RESPECTIVELY. 

• V.X. A.o: EXPT Uq, u e .T w . i 1/2. RESPECTIVELY. 

M: MULTIPLE-TIME-SCALE TURBULENCE MODEL. A : ALGEBRAIC STRESS TURBULENCE MODEL 

(a) Development of the wall- jet flow along the flow direction. 
Figure 4. A wall- jet flow. 
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: COMPUT. RESULTS FOR VELOCITY. TURBULENT KINETIC 

ENERGY, AND REYNOLDS STRESS, *, o, A : EXP'T DATA FOR VELOCITY, TURBULENT KINETIC 
ENERGY, AND REYNOLDS STRESS, M: MULTIPLE-TIME-SCALE TURBULENCE MODEL, 

A: ALGEBRAIC STRESS TURBULENCE MODEL. 

(b) Velocity, turbulent kinetic energy, and Reynolds stress at x = 1.1532 m. 



: kf/kp, : Et'ep, : £t FOR MULTIPLE-TIME-SCALE 

x: DISSIPATION RATE FOR ALGEBRAIC STRESS TURBULENCE MODEL 


(c) Dissipation rate, k^/kp, and £jV e p at x = 1.1532 m. 
Figure 4. (Continued). 
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COMPUT. RESULTS FOR VELOCITY, TURBULENT KINETIC 

ENERGY. AND REYNOLDS STRESS. *, o. A : EXPT DATA FOR VELOCITY. TURBULENT KINETIC 
ENERGY. AND REYNOLDS STRESS, M: MULTIPLE-TIME-SCALE TURBULENCE MODEL. 

A: ALGEBRAIC STRESS TURBULENCE MODEL. 

(d) Velocity, turbulent kinetic energy, and Reynolds stress at x = 1.6882 m. 



: Vkp, : £t/£p, :£t FOR MULTIPLE-TIME-SCALE 

x: DISSIPATION RATE FOR ALGEBRAIC STRESS TURBULENCE MODEL 


(e) Dissipation rate, kj./kp, and e^/Cp at x = 1.6882 m. 
Figure 4. (Concluded). 
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3-4. Wake -Boundary Layer Interaction Flow 

The last example case considered herein is a weakly coupled wake- boundary 
layer interaction flow, the experimental data of which can be found in Tsiolakis, 

Krause, and Muller [18]. The wake was generated by placing a cylinder inside a flat 
plate boundary layer flow. The free stream velocity of the boundary layer flow was 
26.4 m/sec, the diameter of the cylinder was 0.01 m, and the cylinder was separated 
from the wall by a distance of twice the flat plate boundary layer thickness. 

The transverse domain extending from y = 0.002 m in the near wall region up 
to y = 0.12 m in the free stream region was discretized using 45 unequally spaced 
quadratic elements; and the flow direction domain extending from x = 0.2 m (x/d = 20) 
to x = 0.86 m (x/d = 86) was discretized using 825 line steps. Approximately 10 
iterations were required for each line- step. 

The computational results of velocity, turbulent kinetic energy, and the Reynolds 
stress at the two down-stream locations are shown in Figures 5(a) and 5(c). The 
computational results of dissipation rate, k^/kp, and e^/Cp are shown in Figures 5(b) 

and 5(d). Again, all of the turbulence quantities as well as the dissipation rates 
obtained by using the two different turbulence models were almost identical; and the ratio 
of e^/e p showed sharp transition in the viscous super layer. It was found that the 

present computational results compared favorably with experimental data, in general. 



. : COMPUT. RESULTS FOR VELOCITY. TURBULENT KINETIC 

ENERGY, AND REYNOLDS STRESS, x, o, A : EXPT DATA FOR VELOCITY. TURBULENT KINETIC 
ENERGY, AND REYNOLDS STRESS, M: MULTIPLE-TIME-SCALE TURBULENCE MODEL, 

A: ALGEBRAIC STRESS TURBULENCE MODEL 

(a) Velocity, turbulent kinetic energy, and Reynolds stress at x = 0.40 m. 
Figure 5. A wake-boundary layer interaction flow. 
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: kf/kp, : Et^Ep, £t FOR MULTIPLE-TIME-SCALE 

x: DISSIPATION RATE FOR ALGEBRAIC STRESS TURBULENCE MODEL 


(b) Dissipation rate, k t /kp, and e t /ep at x - 0.40 m. 



; COMPUT. RESULTS FOR VELOCITY. TURBULENT KINETIC 

ENERGY, AND REYNOLDS STRESS, x, o, A : EXPT DATA FOR VELOCITY. TURBULENT KINETIC 
ENERGY, AND REYNOLDS STRESS, M: MULTIPLE-TIME-SCALE TURBULENCE MODEL. 

A: ALGEBRAIC STRESS TURBULENCE MODEL 

(c) Velocity, turbulent kinetic energy, and Reynolds stress at x = 0.86 m. 

Figure 5. (Continued). 
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x: DISSIPATION RATE FOR ALGEBRAIC STRESS TURBULENCE MODEL 
(d) Dissipation rate, k t /k p , and e t /e at x = 0.86 m. 

Figure 5. (Concluded). 

IV. CONCLUSIONS AND DISCUSSION 


A set of accurate computational results for complex turbulent boundary layer 
flows, i.e., the wall-jet and wake-boundary layer interaction flows, has been pre- 
sented in this paper. The success may be due to the turbulence models used as well 
as the accurate numerical method which is free of numerical diffusion. 

It was observed that the multiple-time- scale turbulence model could resolve 
details of the complex turbulent boundary layer flows somewhat better than the 
algebraic stress turbulence model. Nevertheless, both of the turbulence models 
yielded equivalently accurate computational results for most of the example cases con- 
sidered, and could resolve the viscous super layer. In this sence, it would be 
appropriate to discuss the relative advantages and disadvantages of the two turbu- 
lence models. 

The advantage of the multiple- time- scale turbulence model lies in its capability 
to model the cascade process of turbulent kinetic energy, its capability to include the 
variation of the effective eddy viscosity coefficient implicitly, and its potential to 
include more experimental observations in the future. 

The disadvantages of the multiple- time -scale turbulence model over the algebraic 
stress turbulence model may seem to be the computational overload to solve two more 
transport equations. But the overload was compensated for by lack of the necessity 
to compute the complicated effective eddy viscosity coefficient (c ) equation used in 
the algebraic stress turbulence model. y 
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APPENDIX I 


THE TURBULENCE MODEL CONSTANTS 


Establishment of the turbulence model constants is based on the assumption that 
the ratios of kjVkp and Dk^/Dk^ will asymptotically approach constant values in simple 

homogeneous turbulent flows such as the decay of grid turbulence [19], a nearly 
homogeneous shear flow [20], and the near wall equilibrium turbulent flows. This 
assumption is justified by the computational results posterily. The model constants 
also need to satisfy the realizability conditions (i.e., both of the turbulent kinetic 
energies, energy transfer rate, and the dissipation rate cannot become negative) 
described below. 


A- 1-1. Free Stream Turbulence 

The velocity gradient and the production of turbulent kinetic energy vanish in 
the free stream region of turbulent boundary layer flows, and the turbulence equa- 
tions can be simplified as: 


Dk 

P _ 

nt 


= - e 


De 

f 

Dt 


- - c 


p3 k 


'P 


P 


Dk t 

Dt 


= £ 


De t 

Dt 


= c 


tl k. 


+ c 


£ p e t 


t2 k 


t 


't3 k. 


(A . 1) 


(A. 2) 


(A. 3) 


(A. 4) 


Let the ratio of the dissipation range turbulent kinetic energy to that of the 
production range be: 


k 
IT 


n = constant 


(A. 5) 


Then, based on the definition of the total turbulent kinetic energy, i.e., 

k = k +k. , it can be obtained that: 

P t’ 
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(A. 6) 


k = 1 

K p (n+1) K 


n 

t ” Tn+TJ 


k 


(A. 7) 


Taking material derivative of equation (A. 5) yields, after rearrangement: 


k t (Dk t /Dt) 
IT = (Dk p /Dt) 


(A. 8) 


Substituting equations (A.l) and (A. 3) into equation (A. 6) yields: 


— = n + 1 
P 


(A. 9) 


It can be seen from equation (A. 7) that the ratio of the dissipation rate to 
transfer rate also becomes a constant in the free stream region. Substituting equa- 
tions (A. 6) and (A. 9) into equations (A.l) and (A. 2), respectively, and dividing the 
results yield: 


Dk - _J_ 

^ " c p3 e t 


(A. 10) 


where the value of c „ is the decay rate of grid turbulence, and is given as, from 

po 

the experimental data [19]: 


c p3 


= 1.8 through 1.92 


(A. 11) 


Substituting equations (A. 7) and (A. 9) into equations (A. 3) and (A. 4) and 
dividing the results yields: 


Dk _ -J_ k_ 

D e. d . £. 

t gt t 


(A. 12) 


where d t = c tl /n/(n+l) + c t2 /n - (n+l)c t3 /n is the decay rate of the grid turbu- 
lence, and hence: 
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(A. 13) 


c tl °J2 _ (n + 1)C t3 = _ 0 

n(n+l) n n p3 


In the free stream region, the eddy viscosity equation can be written as: 


v 


t 



f (n+1) f- 

t 


(A. 14) 


where the constant (n+1) plays the role of the c^ function in the algebraic stress 

turbulence model, and the constant n can be determined from experimental data or 
from the c^ function of the algebraic stress turbulence model. 

A- 1-2. Near Wall Equilibrium Flow 

In the near wall region, where the turbulence is in equilibrium (the production 
rate is approximately equal to the dissipation rate), it can be assumed that: 


Pr = 



e t 


(A. 15) 


As in the standard wall function analysis [13], it is assumed that the convec- 
tive transport of turbulence quantities is negligible compared with the diffusive 
transport, and that both of the turbulent kinetic energies are almost constant in the 
near wall region. Then the turbulence equations can be simplified as: 


, v, 3 k 

^ [(v + 4 )l ^ 1 = ° 


(A. 16) 


j v, 3 c 

ay 1<v + 


<c pi + C P 2 - c p3> 


(A. 17) 


j V » 9 ■« > 

d r/- t \ kt , _ n 

3y ° kt dy 


(A. 18) 


v t 3 

r et 


37 E(v + W 1 ” tS~ (C tl + c t2 " c t3 } 


(A. 19) 


From the mixing length assumption, it can be obtained that: 
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C 3 / 4 k 3/2 c 3/4 (k +k ) 3/2 

e = _H = _H E_J 

<y Ky 


(A. 20) 


v = c — = c — — 

t ye ye. 


vv 


(A. 21) 


where e (e=£p=e t ) is the dissipation rate in the near wall region. Substituting equa- 
tions (A. 20) and (A. 21) into equation (A. 17) yields: 


k 17T7 7 

^ ep y p3 pi p2 


(A. 22) 


In the same way, substituting equations (A. 20) and (A. 21) into equation (A. 19) 


yields : 


1/2 ; 

°£t y (c t3 


(A. 23) 


The ratio of the turbulent kinetic energies, k^/k^, obtained from the energy 

transfer equation and the dissipation equation should be the same in the near wall 
region. Equating equations (A. 22) and (A. 24) yields: 


a (c„-e,-c 0 )=o . ( c . 0 - c , . - c, 0 ) 
ep p3 pi p2' et t3 tl t2 y 


(A. 24) 


Equations (A. 22) and (A. 23) provide two constraint conditions for the turbu- 
lence model constants to satisfy, since neither of the turbulent kinetic energies can 
become negative. 


A- 1-3. Homogeneous Turbulent Shear Flow 

The homogeneous turbulent shear flow experiment due to Harris et al. [20] is 
considered in the following. In the experiment, the turbulent flow with constant mean 
velocity gradient approached asymptotic state in which scales and turbulent kinetic 
energy grow monotonically . The experimental data at the downstream end of the test 

2 3 2 2 

section are given as e. = 3.35 m /sec , k = k + k. = 0.41 m /sec , and Pr = 5.84 

3 2 1 P 1 

m /sec , where the production rate has been estimated using the relationship that 
Pr = -u'v' (9u/3y). For the turbulent flow with constant mean velocity gradient, 
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the production range turbulence equations, equations (3) and (5), can be written as: 


Dk p 

nt - 


Pr 


(A. 25) 


De 
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Dt 


Pr' 


Pre 


+ c 


P 2 K 



(A. 26) 


and the dissipation range turbulence equations can be written as: 


Dk t 

DF 


= e 


(A. 27) 


Dc t 

Dt 


= c 


tl k 


t 


+ c 


Vt 


t2 k 


t 


c t3 k^ 


(A. 28) 


It is assumed that the partitioning of the turbulent kinetic energy spectrum is 
dependent on the ratio of production rate over dissipation rate. Let the ratios of 
production rate to dissipation rate and the energy transfer rate to dissipation rate 
be a and 3, respectively, i.e., 


Pr 

e t 


- a 


(A. 29) 


i = e . (A. 30) 

e t 


Substituting equations (29) through (31) into equations (25) and (26) yields, 
after a few algebra: 


t Dk 

k De t Y 


3( or (3) 


(C pl“ 2 + C p2 ag ' C p3^ 2) 


(A. 32) 


where the constant y can be obtained from the experimental data as shown below. 
The length scale £ is related to the dissipation rate as: 
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(A. 33) 


c k 


3/2 


Taking the derivative of equation (A. 33) yields, after some rearrangement: 


dfc 3 dk de 

Z ^ TT” e 


(A. 34) 


From the experimental data, we have: 


aifx- = °- 6818 


(A. 35) 


Substituting equation (A. 35) into (A. 34) yields: 


k *^ e t 

* =7:-air= °- 82 


(A. 36) 


In the similar way, it can be obtained from equations (27) and (28) that: 


e t dk 
k T 


(a- 6) 


( c ti6 + 


C t3 ) 


(A. 37) 


where the physical implication of y is the same in equation (A. 32). 

The turbulence model constants need to satisfy equations (A. 10), (A. 13), (A. 24), 
(A. 32), and (A. 37) as well as the constraint conditions imbedded in equations (A. 22) 
and (A. 37). 

The turbulence model constants given previously were obtained by assuming that 
the decay rate of grid turbulence (c „) is equal to 1.84, the ratio of k^/k^ in the 

free stream region (n) is equal to 1.5, and the ratios of Pt/e^ and e p/ e ^ f° r the 

homogeneous shear flow [20] are equal to 1.5 and 1.05, respectively. The constant 
was obtained to be 0.29 through computer optimization. Note that the ratio of 

Pr/e^ is equal to 1.74 according to the experimental data given in Harris et al. [20]. 

But the ratio of Pr/e. is equal to 1.74 had been obtained only at the downstream end 

of the test section, and the data would be subjected to experimental uncertainties. 

As more experimental data become available, the turbulence model constants could be 
better determined in the future. 
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APPENDIX II 


CONTROL-VOLUME BASED FINITE DIFFERENCE COMPUTATION 
OF SEPARATED AND SWIRLING TURBULENT FLOWS 


A control volume based finite difference computation of elliptic turbulent flows 
is presented in this Appendix. In turbulent flow computations, numerical uncertainties 
such as grid size, grid spacing, wall function, and numerical diffusion may mask the 
performance of a turbulence model. The details on the control- volume based finite 
difference code (TEACH-T) used herein can be found in Reference 23. The example 
problems considered were a backward- facing flow [22] and a confined swirling jet 
flow [ 23] . 


A- 2-1. A Turbulent Backward-Facing Step Flow 

For the backward- facing step flow case, the upstream boundary has been 
located at five step heights upstream of the expansion corner. Flat profiles for the 
flow direction velocity (U Q ) of 18.5 m/sec, turbulent kinetic energy of 1.027 m^/sec*. 

and the dissipation rate of 780 m^/sec^, have been used as the inlet boundary condi- 
tions. The turbulent kinetic energy and the dissipation rate have been obtained using 

Of o 

the relationships that k = 0.003 U Q 2 and e = c yf k U, respectively, where % 

(= 0.003H) is the length scale [24]. The inlet boundary condition for the ratios of 

k /k and zjz have been assumed to be 0.25 and 1.0, respectively, based on the 
t p t p 

computational result of the fully developed channel flow . The Reynolds number based 
on the inlet flow condition was 45,000. The computed streamline and the turbulent 
kinetic energy contours, for the expansion ratios of 1:2, are shown in Figures A.l 
and A. 2, respectively. The predicted reattachment length was 6.3H, where H is the 
step height, which is in much closer agreement with the experimentally observed value 
of 7(±1)H compared with the 5.2H predicted by using the standard k-e turbulence 
model [ 24] . The computed mean velocity profiles at three cross- stream locations 
(x/H = 5.2, 10.7, and 16) are compared with experimental data in Figure A. 3. The 
velocity profile at x/H = 5.2 indicates that the recirculation zone is smaller than that 
of the experimental data [ 24] . The turbulent kinetic energy (k = k p + kp and the 

shear stress profiles at three cross-stream locations are shown in Figures A. 4 and 
A. 5, respectively. These figures illustrate that the peaks and shapes of these two 
turbulence quantities have been successfully predicted by the present turbulence 
model. It can be found in Reference 24 that the standard k-e turbulence model 
severely under-predicted the reattachment length and the magnitudes of these turbu- 
lence quantities. It can be seen from these computational results that the multiple- 
time- scale turbulence model can yield significantly improved predictions for separated 
turbulent flows than the standard k-e turbulence model. 

The ratios of dissipation rate to energy transfer rate at the same cross- stream 
locations and at x/H = 27 are shown in Figure A. 6. It can be seen that the ratio is 
smaller in the high turbulence region than in the low turbulence region. Thus, the 
turbulent eddy viscosity coefficient is smaller in the high turbulence region than in 
the low turbulence region, see equation (10). At far downstreams, i.e., x/H > 27, 
the flow approached fully developed channel flow state and the ratio (e^/e ) became 
almost uniform across the channel width. 
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The streamline and the turbulent kinetic energy contours obtained by using 
= Oj^ = 1.0 and = 1.3, which are the most frequently used coefficients 

in most of the turbulence models [24], are shown in Figures A. 7 and A. 8, respec- 
tively. These contour lines compare favorably with those in Figures A.l and A. 2, 
except that the reattachment length was slightly shorter than the previous case. 


A 



B 

C 

D 

E 

F 

G 

H 


Figure A-l. Streamline contour for backward -facing step flow. 



Figure A- 2. Turbulent kinetic energy contour for 
backward- facing step flow. 


24 















A- 2- 2. A Confined Swirling Jet 

For confined swirling jet computations, the inlet boundary conditions were 
obtained directly and/or by curve-fitting the experimental data [23] at 0.005 m down- 
stream of the expansion corner in order to reduce uncertainties in the inlet boundary 
condition. The inlet boundary condition for the ratios of k./k and e./e has been 

obtained by using equation (19). A 42 by 30 grid with concentrations of nodes in the 
recirculation regions, in the entrance region, and near the wall was used for the 
calculation. The streamline and the turbulent kinetic energy contours are shown in 
Figures A. 9 and A. 10, respectively; and the development of the axial velocity along 
the center line is shown in Figure A. 11. The predicted strength and size of the 
central recirculation zone compared favorably with the experimental data. The loca- 
tion and shape of the central and corner recirculation zones shown in Figure A. 9 
compared more favorably with experimental data than those obtained by using the 
standard k-e turbulence model and the two other k-e type turbulence models which 
include the swirl-related modifications [24], Comparison of the calculated mean axial 
and swirl velocity profiles with the experimental data are presented in Figures A- 12 
and A- 13, where it is shown that the calculated mean velocity profile mimic the exhibited 
data trend and that these velocities are in good agreement with measurements at most 
of the locations. The predicted axial velocity profiles were skewed toward the wall, 
the trend of which is in good agreement with that of the measured data. At further 
downstream locations, the computed swirl velocity exhibited a slightly premature decay 
toward a forced vortex profile. The discrepancies in these predictions can be partially 
attributed to the wall functions used. The free vortex swirler used in the experimental 
study [23] generated a very unstable flow field. The equilibrium assumption under- 
lying in the wall functions may be doubtful in this situation. Nevertheless, the suit- 
ability of the present turbulence model for complex turbulent flows has been partly 
demonstrated through these computational results. 

2 

The r.m.s. values of the fluctuating axial velocity ( u 1 ) at various downstream 
locations are compared with experimental data in Figure A- 14. The computational 
results showed similar trend with that of experimental data, but the magnitude of the 
r.m.s. values of the fluctuating axial velocity was found to be under-predicted. 

These under-predictions may be due to the isotropic turbulence assumption used in 
estimating the r.m.s. value of the fluctuating axial velocity from the computed turbu- 
lent kinetic energy. 

The ratios of dissipation rate to energy transfer rate at the same axial loca- 
tions are shown in Figure A- 15. It can be seen that the ratio is smaller in the high 
turbulence region than in the low turbulence region, which was the same for all the 
flow cases considered herein. At far downstreams, the flow approached fully developed 
pipe flow state and the ratio (e^/e^) became almost uniform in the radial direction. 

The streamline and the turbulent kinetic energy contours obtained by using 
a, = a. . = 1.0 and a = a . = 1.3 are shown in Figures A- 16 and A- 17. Only very 

Kp Kt ep Cl 

insignificant changes were detected in mean features of the flow such as the stream- 
line contour and mean velocity profiles. Slight difference was found in the turbulent 
kinetic energy contours. The computational results obtained by using a. = =1.0 and 

a = a . = 1.3 were found to be more diffusive than those obtained by using the 

sp ct 

other set of coefficients, however, the differences were insignificant in general. 
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Figure A- 12. Axial velocity profiles of the confined swirling jet 
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Figure A- 13. Swirl velocity profiles of the confined swirling jet 
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TABLE A- 1 . CONTOUR LABELS FOR BACKWARD-FACING STEP FLOW 


(a) Streamline (^*) 


Label 

* 

Label 


Label 


A 

0.0 

E 

0.4 

I 

0.6081 

B 

0.1 

F 

0.5 

J 

0.685 

C 

0.2 

G 

0.6 

K 

0.69 

D 

0.3 

H 

0.6767 

L 

0.695 


: The streamline data are based on the normalized velocity. 


(b) Turbulent Kinetic Energy (k**) 


Label k 

Label k 

Label k 

A 0.0394 

B 0.0365 

C 0.0292 

D 0.0219 

E 0,0146 

F 0.0073 


** : The turbulent kinetic energy has been normalized by Uq^ p 
where Uq ( — 18.5 m/sec) is the inlet velocity. 


TABLE A-2. CONTOUR LABELS FOR CONFINED SWIRLING JET 

(a) Streamline (V>*) 


Label 

i> 

Label 

■ * 

Label 

v> 

A 

-0.25E-4 

F 

0.80E-4 

K 

0.60E-3 

B 

-0.75E-5 

G 

0.12E-3 

L 

0.62E-3 

C 

-1 . 0E-7 

H 

0.24E-3 

M 

0.65E-3 

D 

0 . 2E-4 

I 

0.36E-3 

L 

0.69E-3 

E 

0.4E-4 

J 

0.48E-3 




* : The streamline data are based on the physical velocity. 


(b) Turbulent Kinetic Energy (k, m^/sec^) 


Label k 

Label k 

Label k 

A 0.15 

B 0.13 

C 0.11 

D 0.08 

E 0.05 

F 0.03 

E 0.01 

F 0.008 
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